1/1 


AD-A172  989 
UNCLASSIFIED 


MINIMIZING  THE  REFLECTION  OF  ELECTROMAGNETIC  NAVES  AV 
SURFACE  INPEDANCE(U)  WISCONSIN  UNIV-HADISON  MATHEMATICS 
RESEARCH  CENTER  T  J  BRIDGES  ET  AL.  JUL  86  HRC-TSR-2942 
DAAG29-88-C-8841  F/G  28/14 


MRC  Technical  Summary  Report  #2942 

MINIMIZING  THE  REFLECTION  OF 
ELECTROMAGNETIC  WAVES  BY 
SURFACE  IMPEDANCE 

T.  J.  Bridges,  G.  Chen  and  G.  Crosta 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 

July  1986 

(Received  June  20,  1986) 


DTIC 

electei 

OCT  8  1986 


OTIC  FILE  COP* 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

U.  S.  Array  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation 
Washington,  DC  20550 


86  10 


UNIVERSITY  OP  WI SCONS IN-MADISON 
MATHEMATICS  RESEARCH  CENTER 


MINIMIZING  THE  REFLECTION  OF  ELECTROMAGNETIC  WAVES 
BY  SURFACE  IMPEDANCE 

1  *  _  ?  **  1  *** 

T.  J.  Bridges'  ,  G.  Chen  '  and  G.  Crosta  ' 

Technical  Summary  Report  #2942 
July  1986 

ABSTRACT 


t 


In  an  empty  halfspace  a  point  source  emits  electromagnetic  waves  of  fixed 
frequency  and  arbitrary  polarization.  The  fields  reflected  by  an  imperfectly 
conducting  plane ,  characterized  by  a  constant,  isotropic  surface  impedance, 
are  determined  by  means  of  the  geometrical  optics  approximation.  An 
optimization  method  is  used  to  compute  the  surface  impedance,  which  minimizes 
a  given  function  of  the  e.m.  fields  (e.g.  the  average  energy  density)  at  a 
given  observation  point.  The  properties  of  the  functions  to  be  minimized  are 
studied  and  a  set  of  numerical  results  is  presented  and  discussed. 
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MINIMIZING  THE  REFLECTION  OP  ELECTROMAGNETIC  WAVES  BY  SURFACE  IMPEDANCE 


T.  J.  Bridges1'*,  G.  Chen2'**  and  G.  Croat*3'*** 

I.  INTRODUCTION 

1.1.  Problem  ititwent 

We  consider  an  optimization  problem  related  to  the  reflection  of  electromagnetic 
waves.  A  point  source  of  monochromatic  waves,  characterized  by  a  dipole  vector,  is 
located  in  an  empty  halfspace  and  illuminate*  an  unbounded  plane  surface,  having  constant 
surface  impedance.  If  the  latter  satisfies  some  properties  to  be  defined  below,  the 
incoming  wave  energy  is  partly  absorbed  and  partly  reflected.  By  means  of  the  geometrical 
optics  approximation  the  reflected  fields  are  evaluated.  Since  they  are  functions  of 
position  as  well  as  other  parameters,  we  state  the  following  problems: 

a)  given  an  observation  point  in  the  half space,  which  need  not  coincide  with  the 
source,  determine  e.g.  how  the  e.m.  energy  density,  or  the  real  part  of  the 
modulus  of  the  Poynting  vector  at  that  point,  depend  on  the  surface  impedance  as 
its  real  and  imaginary  part  are  varied  (direct  problem): 
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b)  given  an  observation  point  as  above,  find  the  value  of  the  surface  impedance 


within  a  prescribed  set  which  minimises  either  of  the  above  mentioned  field 


functions  (minimisation  problem); 


c)  after  finding  the  optimal  impedance  value  as  of  problem  b) ,  determine  how  either 


of  the  above  field  functions  depends  on  position  as  the  observation  point  is 


moved  (position  sensitivity  problem). 


Problem  a)  consists  of  making  the  reflecting  plane  almost  'invisible*  to  a  particular 


observer.  This  study  is  our  first  approach  towards  the  widely  publicised  "stealth" 


problem,  which  aims  at  minimizing  the  radar  cross  section  of  some  obstacles.  The  model  we 


adopt  is  simplified,  however,  with  respect  to  any  practical  situation:  nevertheless  we 


believe  our  methods  and  results  are  helpful  in  outlining  some  features  of  the  general 


problem. 


In  the  following  we  introduce  the  background  material.  In  Section  II  we  calculate 


the  first  two  terms  in  the  geometrical  optics  expansion  of  the  reflected  magnetic  and 


electric  fields.  In  Section  III  we  state  the  optimization  problems  for  some  objective 


functions,  which  depend  on  the  complex  surface  impedance  through  the  reflected  fields. 


Mnimlzatlon  of  said  functions  is  carried  out  by  a  steepest  descent  method,  incorporated 


into  a  recently  developed  computer  code.  Several  results  are  shown  and  discussed. 


1.2.  The  impedance  boundary  condition 


Let  D  be  a  halfspacei  x  -  (x1(x2,Xj),  where  x3  >  d,  denotes  a  point  in  D  and 
3d  the  x3  -  d  plane.  Then  5  will  stand  for  D  U  3d.  Wie  electromagnetic  fields 


(H,B)  considered  herewith  are  complex  quantities  having  a  suppressed  time  dependence 


exp(-lut).  Both  H  and  K  satisfy  homogeneous  Maxwell's  equations  everywhere  in  D, 


except  at  Xg,  where  the  source  of  a  monochromatic  spherical  wave  is  located  (see  Section 


1.5).  We  assume  that  the  x3  «  d  plane  is  an  imperfect  conductor.  In  order  to  account 


for  the  nonzero  resistance  of  the  latter  and  for  the  related  inductive  effects,  the 


following  impedance  boundary  condition  is  specified: 
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-m**«Z(i»*nxH),  at  3d 
where  a  ia  the  outward  unit  nornal  vector  and  Z  ia  a  complex  valued  3*3  matrix, 
named  the  eurface  Impedance.  In  the  moat  general  case  Its  entries  may  depend  on  i  e  >D. 
It  can  be  shown  (Jones  [1],  pp.  S11  ff)  that,  if: 

Z  +  Z*  >  0  ,  (1.2) 

where  Z*  ia  the  hermitian  conjugate  of  Z,  i.e.  if  Z  ♦  Z*  is  positive  definite,  then 
(1.1)  plus  some  adequate  radiation  conditions  are  sufficient  for  Maxwell's  equations  to 
have  a  unique  solution  in  D. 

If  Z  is  replaced  by  the  complex  scalar  i  i>  p  «  it,  (1.2)  becomes: 

p  -  Re  z  >  0  ,  (1.3) 

which  is  easily  seen  to  yield  a  power  dissipative  boundary  condition. 

Before  considering  our  problems  in  detail,  we  must  conusant  on  (1.1).  It  can  e.g.  be 
obtained  under  the  hypothesis  that  fields  penetrating  from  the  surface  into  the  material 
fall  off  as  (Collin  [2],  pp.  16-17)  e.g.: 

K(m)  -  Z(x1,x2,d)  exp[-(d  -  x3)/«8),  x3  <  d  (1.4) 

where  6S,  the  "skin  depth",  is  related  to  conductivity  a,  angular  frequency  u  and 
vacuum  magnetic  permeability  lig  by: 

Sg  :■  /2/ii>UqO  .  (1.5) 

If  the  material  is  isotropic,  it  can  be  shown  that  the  corresponding  surface  impedance 


*  -  (1  +  i)/crSs 

Moreover,  for  (1.5)  to  hold  it  is  required  that: 


«.  <<  * 


w/Vo 


(1.6)  and  (1.7)  respectively  shall  be  used  as  a  physical  realizability  constraint  in 
the  minimization  problem  and  as  an  "a  posteriori"  inequality  for  the  validity  range  of  our 
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.  The  geometrical  optics  solution  method 


According  to  geometrical  optics  and  diffraction  theory  (Luneburg  [3],  Keller  [4], 


Klein  and  Kay  [S] )  the  asymptotic  solution  to  a  diffraction  problem  may  be  expressed  as 


the  sum  of  the  incident  and  the  scattered  fields i 


<B,B)  -  (B1,*1)  +  (B®, B®)  , 


where  the  scattered  field  (BS«BS)  is  further  decomposed  Into  the  sum  of  the  reflected 


and  diffracted  fields: 


(B8,*3)  -  <Br,Br)  ♦  (B0,!0)  . 


All  three  fields  satisfy  Maxwell's  equations*  In  our  analysis  we  assume  that  there  is  no 


diffraction,  l.e.  (BD,BD)  «  (0,0).  We  then  postulate  that  the  reflected  field  has  a 


uniformly  valid  asymptotic  expansion  of  the  form: 


(BR(x),BR(x))  ~  eikS(x)  l 


(B(n)(x),B(n)(x)) 


(1.10) 


n-0  (ik> 


where  k  :«  u/c  and  c  :«  1/(e0UQ)  ^2.  (BR,BR)  is  an  outgoing  wave  satisfying  the 


homogeneous  Maxwell's  equations  and  radiation  conditions. 


It  is  easily  shown  that  S(*),  B*n\  B*n)  satisfy  respectively: 


|Vs|2  -  1  (eikonal  equation) 


(1.11) 


and,  along  a  ray: 


2<j^B(n)  +  (V2S)B*n*  -  _v2*<n_1) 


(1.12) 


2dyH(n)  +  (V2S)B(n>  -  -V2B<n-1), 


n  -  0,1,2,... 


(1.13) 


with  H' 


0 ,  where : 


d^xdl)  =  VS(x)  •  Vf  (x(Y ) ) 


(1.14) 


is  the  differentiation  along  a  ray. 


The  boundary  condition  (1.1)  becomes: 


B^  —  B(J 


H*)  -  n  *  B1  -  B(i 


H1)  at  3D 


(1.15) 


In  some  special  cases,  e.g.  with  TK  and  TM  waves.  Maxwell’s  equations  reduce  to  6 


scalar  wave  equations  with  uncoupled  impedance  BCs  (Senior  [6]).  This  was  the  case  with 
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an  earlier  worX  of  ours  (Chan  t  Bridges  [7]),  where  we  considered  a  problem  for  acoustic 
waves.  For  tha  praaant  problan,  no  such  raductlon  la  possible,  hanca  tha  calculations 
will  ba  laora  complicated. 

1.4.  A  preliminary  example t  a  plana  wavs  lncldant  on  an  Impedance  plana 

To  illustrata  tha  basic  rola  of  an  impedance  boundary  condition,  tha  following  simple 
example  is  considered.  An  E-polarised  wave  R1  i-  (0,0, Rj)  is  incident  on  the  (x1,x3>- 


plane ,  where : 


T  -iMx.cosa+x0sina ) 

Rj(x)  -  e  1  2  ,  (a  j* 


(1.16) 


a 


I 


Figure  1.  A  plane  wave  lncldant  on  an  impedance  plana. 


The  Xj  "  0  plane  is  assumed  to  have  the  following  boundary  impedance > 


‘s  0  o" 

0  s  0  , 

.0  0  t. 


(1.17) 


A  straightforward  calculation  based  on  Sections  1.1  snd  1.2  gives  the  following  exact 
solution  for  tha  reflected  field i 


„  „  R  V  ”  **  *in°  -ih(x.coso-x2sinB ) 

*l(x)  •  *2<*>  “  “  VJ V  • 


(1.18) 
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u  Qu  >  Ice  Etna 


h5(x)  -  0  . 

From  the  last  expressions,  the  special  value i 


x  «  u0c/(aina) 


<1.19) 


<1.20) 


implies  that  the  reflected  fields  vanish  everywhere  in  Rj.  This  example  helps  us  to 
point  out  some  facts,  which  hold  in  a  sore  general  context i  (a)  the  optimal  Impedance, 
which  minimizes  the  reflected  field  depends  upon  the  angle  of  incidence  of  the  incoming 
wave  and  (generally)  on  its  frequency,  (b)  if  we  l3t  a  «  i/2  in  (1.20),  we  obtain 
z  -  (Uo/eQ)'^2*  which  does  not  comply  with  (1.6),  and  (c)  if  we  assume 
Ra  t  •  (Ug/eg)1/2  :«  ig  ( the  "free  space"  impedance.  Eg  ~  366.9  Ohm),  from  (1.5)  we 
get  a  »  eu/2  and  &a/X  *  1/*,  i.e.  the  wavelength  and  the  penetration  depth  are 
comparable,  contrary  to  the  constraint  (1.7). 


1.5.  A  spherical  electromagnetic  wave  impinging  on  an  Impedance  plane 

Let  an  electric  dipole  antenna  with  moment  p  ”  (p^PjrPj)  be  located  at  the  source 
point  >£  «  (0,0, 2d),  d  >  0.  An  impedance  plane  is  located  at  x3  •  d.  The  reference 
frame  and  the  field  layout  are  shown  in  Figure  2.  On  the  x3  -  d  plane  the  impedance 


boundary  condition: 


-n  »  K  -  z(n  *  n  *  H) 


(1.22) 


subject  to  (1.3)  is  satisfied.  This  is  the  example  we  shall  consider  with  greater 


detail.  In  order  to  solve  any  of  the  minimization  problems  listed  in  Section  1.1,  we  need 
the  reflected  fields  at  the  observation  point  e  D.  These  fields  will  be  evaluated  as 
explained  in  Section  1.3. 


-*•  '  V.V.VV 


*  ■/  /  *„*  %*  %■  •.*  ».  %•  •/  V  »„■  *  *  •  *  s  • 

A.'  C  *l'  v'  •.  "«.*  V*  v'l",- 


XX.  GEOMETRICAL  OPTICS  APPROXIMATIONS 


XI. 1.  The  potential  and  the  Incident  fields 
The  vector  potential  A  is  defined  by: 

V  x  A  -  U0B  , 

(2.1) 

7  •  A  ■  iuiigSo#  (Lorentz  gauge  condition)  , 
where  the  scalar  potential  *  satisfies: 

74  «  -B  +  Lit  A  .  (2.2) 

Then  A  and  *  satisfy  a  Helmholtz  equation,  in  particular  A1,  associated  to 
(•^H1),  satisfies: 


f  V2Ar  ♦  *2Ar 


"4* pi  (X  -  *g) 


(2.3) 


i  radiation  conditions  • 

The  solution  to  (2.3)  (fundamental  solution  or  free  space  Green’s  function)  reads: 

A1  -  peikrVr*  ,  (2.4) 

where  r*  :«  |x*|  »  |x  -  x8|  -  (x2  *  x|  +  (*3  -  2d ) 2 ) 1  / 2  and  p  «  -iu>u  Qp/4*  .  We  are  thus 
modelling  a  dipole  antenna  emitting  spherical  waves.  The  associated  incident  fields  at 
any  point  x  e  D,  x  f  xa,  are: 


H1  -  (t/u0)7  «!«(**  pH«81tr*/w0r*)(1  -  1/ikr*)  (2.5) 

B1  -  (i/eQu)7  x  h1  ,  (2.6) 

where  x*  :=»  (x  -  x8)/|x  -  x8|  “  x*/r*.  (See  e.g.  JacKson,  [8],  $9.2). 


II. 2.  The  boundary  condition 

Since  the  reflected  wave  is  also  spherical,  then  by  Snell's  law  it  diverges  from  a 
virtual  point  source,  which,  as  Figure  2  shows,  is  located  at  the  origin.  Two  geometries, 
spherical  and  planar,  are  involved  in  the  boundary  condition.  Since  n  •  (0,0,-l), 

(1.15)  and  (1.22)  yield  two  scalar  equations: 

*"l  -  *2  “  -*H1  +  «2  *“  b1 
zHj  -  B^  «  -zldj  -  b2  , 

which  are  not  sufficient  to  find  (■**,  HR) .  They  must  be  completed  by  an  independent 


(2.7) 
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scalar  aquation,  l.a. 


■*-<>. 


(2.8) 


Sinca  b1  and  b2  In  (2.7)  ara  known  functions,  obtalnad  from  (2.5)  and  (2.6),  and  slnca 
«*  can  ba  written  in  tanas  of  ■*,  (2.7)  bacaaast 


at  *3  -  d,  VXf,x2 


lksh,  +  *q(^ 3h1  "  31h35  *  bi 
lkzh2  ■  ro(^2h3  "  ®  3h2^  *  b2  1 

where t 

3^  denotes  the  partial  derivative  w.r.  to  x^,  1  -  1,2,3; 
h^  denotes  the  reflected  magnetic  field  components; 

Cg  has  been  defined  in  Section  1.4. 

II. 3.  Geometrical  optics  expansion 

Expansion  (1.10)  for  the  reflected  magnetic  field  becomes: 


(2.9) 


V 

m 

V 

h2 

-  ik  eikt  l  (ik)“n 

vn 

n-0 

.*3. 

-wn. 

r  >  d  , 


(2.10) 


since  S(x)  -  r,  the  coordinate  along  a  ray  being  r.  For  convenience  we  have  factored 
lk  out  of  the  sum.  By  (2.10)  a  system  of  partial  differential  equations  (2.8,  2.9)  is 
transformed  into  an  algebraic  system,  which  is  solved  by  the  Identity  principle.  The 
expressions  of  b^  and  h2  are  polynomials  containing  powers  of  (ik).  Then  the  sets  of 
equations  for  (  -n.,rn«wn)  will  be  obtained  by  equating  the  coefficients  of  given  powers 
of  (ik).  Let  us  begin  with  n  -  0:  by  introducing  polar  coordinates  (see  caption  of 
Figure  2  for  the  notations)  the  following  system  is  obtained: 

(SgCosO  +  s)Ug  -  sosin0cosewg  »  b01 


(atgCosB  +  s)Vg  -  CgSinSsinswg  ’  b02' 
sinBcosfUg  +  sin8sin»Vg  +  cos6wg  -  0 


d,  Vx.,x 


1  '2 


(2.11) 


where: 


-9- 


b01  (1/Wgrp)(XgCOs8  -  r)  (p3sln9»lrup  +  pjCoaB  )  +  ZgSin^ cos»  (PjCoas  -  plains  )i 

b02  <  Vu0rp)(z  “  *0cos8)(p3sin8cos*  +  Pico8®>  +  tgsin^sin*  (p2cosf  -  p^sir*  )  i 

rp  s-  (x2  +  xij  +  d2)1^2  -  d/cos9  . 

Prom  (2.11)  we  get  the  Oth  order  component  of  the  reflected  magnetic  field  at  x3  -  d, 
Vx1 ,x2: 

U0  -  (OQ  -  SPX)/y0rpV»  V0  -  -(UR  +  TPX)/u0r0V,  wQ  -  YP  sine/u0rp  (2.12) 

where i 

P  p^slnp  -  p2cos«  Q  p3sin8sine  ♦  pjCos8 

R  :«  p3sin8cose  *■  p^coeS  S  :»  ZgSin28coa<p 

T  :«  ZgSln  Seine  U  s-  ZgcosS  -  z 

V  :«■  ZgCoa8  +  z  X  s-  2z  cos8/(Zg  +  z  co»8  ) 

Y  i-  (zQ  -  z  coe8)/(z0  +  z  coed)  .  (2.13) 

The  result  stated  by  (2.12)  can  be  extended  Vx3  >  d  by  Integrating  the  transport 

equation  (1.13).  This  is  a  straightforward  step  for  n  -  0  only,  since  then  the  right- 
hand-side  is  0.  The  required  field  components  are  obtained  by  replacing  rp  by  r  in 
(2.12),  where  r  (x2  ♦  x|  ♦  x2)1^2,  Xj  >  d. 

The  Oth-order  reflected  electric  field  components  are  derived  from  (2.12)  and  the 
counterpart  of  (2.6),  by  applying  the  identity  principle.  In  analogy  with  (2.10)  they  are 
denoted  by  (5g.hg.Cg)  and  readi 

50  -  Zg(v0cos8  -  w0sin8slne) 

nQ  -  -Zg(UgCos8  -  WpSinScoeg),  Vx  r  D  (2.14) 

C0  -  ZgSln8(UgSine  -  Vgcos?)  . 

In  order  to  obtain  higher  order  components,  we  may  follow  the  already  mentioned 
procedure  of  applying  the  identity  principle  to  the  algebraic  system  and  solving  the 
corresponding  transport  equations.  This  leads  to  very  complex  calculations.  An 
alternative  approach  is  suggested  by  the  symmetry  properties  of  the  generalized  spherical 
waves  appearing  in  (2.10).  By  applying  a  result  due  to  Keller,  Lewis,  Seckler  [9),  it  can 


be  easily  shown  that 


<1/2r)B 


(2.15) 


where  the  Beltrami  operator  B( • )  is  defined  by: 


B(  • )  (sin8)"1(aflsinea9(*  )  +■  (1/sin8)32  (•  )) 


(2.16) 


In  oar  case  it  can  be  shown  that  none  of  the  quantities  dealt  with  are  singular  at  6-0. 


The  action  of  B(  • )  on  the  components  given  by  (2.12)  yields: 

bu0  -  a|9u0  -  { (zqQ  +  p(s'9x  +  s* Xg) loose  -  2x(s"p?  +  s-p)  +  oq) 
•  (l/u0rv)  +  z0(UQ  -  SPX)/w0rV2  , 

Bw„  -  a2w„  ♦  ( ( KoP  +  YPa )cos8  -  YP  sin^l/u-r  , 


(2.17) 


(2.18) 


where  the  subscripts  8  and  e  denote  the  partial  derivative  with  respect  to  said 
variables,  the  prime  (')  denotes  division  by  sin6  and  the  double  prime  (")  division  by 
sin28.  similarly  BvQ  is  obtained  from  (2.17)  by  orderly  replacing  (Uq,Q,S)  by 
(v0,-R,T)  and  their  partial  derivatives. 

Finally,  the  1st  order  components  of  the  reflected  electric  field  are  written  as 
functions  of  the  already  calculated  0th  and  1st  order  reflected  magnetic  field 
components  and  of  their  partial  derivatives  with  respect  to  polar  coordinates: 


Ao 

it(r,e,e)  *  / r-  (v1 


cos8  -  w.sinSsine  +  cos8 


3v  .  a  3v 
_ 0  sin  8  0_ 

3r  r  38 


3w„  ,  3w„  3  w. 

...  0  cos6  sine  0  cose  0. 

-  smesine  ^ - ; - ^  3— ) 

f^O  2u0  gin9  2uq 

)^(r,8,e)  -  /  —  (-u.,co88  +  w^sinScose  -  cos8 


+  sinScose  + 


3w_  .  3w_  .  3w 

0  cosB  cose  0  sine  0 


38  r  sin8  5e"' 


R 

^(r.e.e)  -  »  (ul 


sinSsine  -  v. sinScose  +  sin8sine 


3  u.  .  ,  3u. 

0  co so  sine  0 


cose  3u0  .  3v0  co s8  cose  3v0  .  sine  3v0, 

+  rsse  TT  '  Bln0co,?  r - W  +  TTuS-  sr*1 


(2.19) 


It  may  be  Interesting  at  this  point  to  make  a  comparison  between  the  asymptotic 
expansion  we  have  chosen  and  other  approximations,  known  as  the  paraxial  and  Fresnel’s 
approximations,  also  frequently  used  in  electromagnetics.  In  our  case  we  have  started 
with  the  fundamental  solution  (2.4)  and  approximated  the  reflected  fields  by  the  Oth 
and  1st  order  terms  in  the  expansion  (2.10),  which  is  acceptable  whenever  kr  >>  1.  0 
the  other  hand  the  paraxial  and  Fresnel's  approximations,  the  equivalence  of  which  has 
been  discussed  by  Crosta  (10]  for  a  scalar  case,  apply  if  the  fundamental  solution  is 
known  in  closed  form  and  affects  phase  factors  appearing  in  the  complex  exponential 
propagator.  Closed  form  fundamental  solutions  are  easily  obtained  when  boundary 
conditions  are  of  the  Dlrichlet  or  Neumann  type,  and  yield  the  well-known  Rayleigh- 
Sommerfeld  propagators.  *nie  solution  corresponding  to  impedance  BCs  for  the  scalar  case 
has  been  given  by  Malyughinets  [11].  In  this  paper  the  implementation  of  Malyughinets' s 
results  will  not  be  considered. 


III.  FINDING  THE  OPTIMAL  SURFACE  IMPEDANCE  VALUE 


III.l.  Problem  statement 

Let  us  assume  that  the  surface  impedance  can  take  on  values  belonging  to  either  of 
the  following  sets: 

zad  !"  {z  4  C  I  0  <  p,  <  Re  z  <  P2,  \t\  <  b}  (3.1) 

and 

Zad  =■  {*  t  I  to  «  *  I*  *)  .  (3-2) 

hereinafter  named  "the  admissible  sets"  and  denoted  by  Za(j,  whenever  no  ambiguity 
arises.  Equation  (3.1)  contains  in  particular  the  positivity  constraint  for  Re  z»  (3.2) 
represents  the  physical  realizability  constraint,  discussed  in  $1.2. 

Let  x  denote  the  index  function  of  Zad,  which  is  0  if  z  «  za<J  and  4» 


elsewhere:  then  we  define  the  objective  functions  Ja,  m  «  1,...,4: 

J,(Z)  -j  (f gl^(x)  •  «,,(*)  +  V0e£(x)  *  HR(*>)  xl*l  —  Wg  +  X  1*1  ,  (3.3) 
J2(z)  Jv  W,dv  +  x 1*1  :-«2+xl*|  .  (3-4) 
J3(*)  «-  ^  |V*>  *  *&*>  I2  +  X  1*1  *-  " 3  +  xl*l  .  (3.5) 
J4(*)  :-  j  Jr  Re [But*)  *  *$,(x)]  •  wff  ♦  X 1*1  «-  *4  ♦  x 1*1  •  (3.6) 


Each  of  the  Jm's  is  the  sum  of  a  "physical"  term,  «*m,  and  of  an  "economical"  one, 

X | *) l  the  purpose  of  the  latter  will  be  clarified  at  the  end  of  this  subsection.  The 
term  W.,  is  the  time  averaged  reflected  energy  density  at  *  “  (r,9,e)  e  D,  while  Wj 
results  from  integrating  M1  over  a  bounded  domain  V  QC  D»  W3  is  the  time  averaged 
squared  modulus  of  the  reflected  Poynting  vector  8g(  x) 1  M4  is  the  flux  of  Re(Sg) 
through  a  given  bounded  surface  1,  s.t.  I  n  3d  ■ 

The  Jffl's  depend  on  z  through  the  fields. 

We  now  state  the  minimization  problem: 

find  the  surface  Impedance  z  t  Zad,  which  minimizes  one  of  the  J^’s,  i*s>: 

J_(z)  “  min  J_(z):  m  given  .  (3.7) 

z.d 


This  problem  will  be  solved  by  an  iterative  procedure,  which  consists  of  the  steps  also 
shown  by  Figure  3: 


1  -  set  k  -  0; 

.k 


2  -  given  1  e  zad'  determine  the  reflected  fields  to  the  desired  order  of  approximation 

( one  ,  in  our  case ) i 

3  -  compute  wjj,  i.e.  jJ|,  since  rk  e  2 sd; 


4  -  from  J*  and  its  derivative  w.r.  to  *,  jJ^,  defined  in  till. 2,  determine 


V+1  ,  7  ¥  Tk+1  .  Jf. 

z  €  Zad  8,t*  Jm  <  Jm» 

5  -  repeat  the  whole  procedure  from  step  n.2  until  some  stopping  condition  is  met. 

In  detail,  step  n.2  requires  the  solution  of  a  control  problem,  described  by  the  map: 


T  :  2 


'ad 


->  (C°*(D))3  *  (c"(D)  )3  :-  X 


(3.8) 


* - >  <V«R> 

where  T  is  a  c“( zad)  dif feomorphism  and  *  is  the  input. 

Step  n.3  is  related  to  observation  of  the  system,  since  Wm(* )  is  an  observation  map 
and  Wffl( t)  is  the  output. 

Step  n.4  consists  of  updating  the  input  according  to  the  observed  quantity;  in  other 
words  information  is  feedback  from  output  to  input. 

Among  the  stopping  conditions  appearing  in  step  n.5  the  following  are  considered: 

-  an  upper  bound  on  the  number  of  iterations,  k 

-  a  lower  bound,  Jm(nd.n)>  Cor  the  objective  function 

-  a  lower  bound,  lJmz(mln)l'  for  tt,e  ra°dulus  the  derivative  (see  §111.2). 

Without  going  into  further  details,  we  point  out  that  both  analytiu^l  and  physical 
considerations  play  a  role  in  determining  Jn(min)*  us  assume  that  the  minimum 

detectable  energy  in  a  given  V  is  known:  this  value  may  then  be  assigned  to 

From  an  analytical  viewpoint,  can  be  related  to  the  modelling  error,  which 

affects  all  field  functions,  since  (2.10)  and  its  counterpart  for  the  electric  field  are 


replaced  by  the  sum  of  the  first  two  terms.  With  reference  to  Wj,  is  chosen  as 


the  best  upper  bound  of  the  energy  density  due  to  the  neglected  terms  (n  >  2  in  (2.10)). 


We  now  consider  the  well-posedness  of  the  minimization  problem. 

Existence  of  at  least  one  z  ig  insured  by  the  economical  tens  defined  above.  Moreover, 

z  continuously  depends  on  the  problem  parameters,  e.g.  the  observation  point  coordinates 

because  T  of  (3.8)  is  smooth  and  all  of  the  J_'s  are  continuous  functions  of  z  in 

m 


Uniqueness  of  z  cannot  be  proved  in  a  straightforward  way.  Other  properties  of  z 
appear  when  the  first  order  derivatives  of  Jm's  are  studied. 


III. 2.  The  derivatives  of  w.r.  to  z 

We  notice  that  all  of  the  Jm's  are  continuous  real  valued  functions  of  a  complex 
quantity,  z,  where  a  field  vector  and  its  complex  conjugate  always  appear  together. 

Since  B(  • )  depends  on  z  alone  and  ■*(*)  depends  on  z*,  then  it  makes  sense  to 

consider  both  z  and  z*  as  independent  variables. 

Let  us  denote  by  ^  the  derivative  of  B  w.r.  to  z :  then  the  first  variation 

of  W^ ( • )  reads: 

dw1  “  \  I€o(**  *  Ezdz  +  *z*  *  *  ♦  M0(«  •  +  h£.  •  H  dz*)]  - 

-  Re(J1z)dp  -  Im(Jl!S)dT  ,  (3. 


J1z  €0I 


*z  U0H  '  "z 


Similar  expressions  are  obtained  for  the  other  physical  terms:  in  particular 

J3z  -  S  •  (H*  *  Hj)  -  S*  •  (■*  »  V  ,  (3. 

J4z  -  <■*  x  H*  +  B*  x  hz)  •  n  ,  (3. 

where  each  of  the  J^'s  is  a  complex  quantity,  hereinafter  called  the  "gradient"  of 
Jm  w.r.  to  z.  The  differential  properties  of  Jm  lead  to  the  following  result: 

let  z  c  2a^,  then  Z  -  0  is  a  locus  of  stationary  points  of  Jm  (for  any  m) , 
which  are  either  relative  maxima  or  minima. 

The  proof  is  straightforward:  from  the  definitions  (3.3)  to  (3.6)  we  get 


J(z)  »  J(z  ),  where  we  have  dropped  the  subscript  m  for  simplicity.  For  a  given  P 


•way  f row  the  boundary  of  z“d 


A 


away  frow  tha  boundary  of  z“d  we  have  J(p>t)  •  J<pi-t)i  aince  J(»  r  )  la  at  leaat  of 
class  C1  w.r.  to  t,  then  dJ(pjt)  -  -dJ(Pi-T):  In  pertlcular  4J(P»0)  ■  0. 

Second  order  analysis  la  needed  in  general  to  tell  whether  at  (Pf0)J(p<0)  attains  a 
local  maximum  or  minimum.  If  we  deal  e.g.  with  J^*)  and  consider  the  analogy  between 
the  reflection  of  spherical  waves  and  the  preliminary  example  of  |l.4,  we  could  easily 
•how,  without  performing  any  computation,  that  real  values  of  a  must  yield  local  minima 
of  J^.  this  intuitive  conclusion  is  supported  by  a  study  of  the  Hessian  matrix  of  J1# 
denoted  by  V2Jj,  and  defined  byi 


V2J1  s« 


tt  “m4 


Jz*z  ‘Tz*z» 


(3.13) 


where  Jzz  stands  for  the  2nd  order  derivative  of  w.r.  to  z.  The  2nd  order 

variation  dJ 1  of  then  readst 


d2J,  -  (ds  dz*) 


alternatively,  it  can  be  expressed  in  terms  of  ( dp ,  dr  )  : 


d2J,  -  (dp  dt)  •  1^ra  •  V2J,  •  *  *^<5P^ 


(3.14) 


(3. IS) 


where  *  t» 


and  T^r-  denotes  the  transpose  of  T.  Let  us  now  define 


Ttr*  •  V2j.,  •  Tf  then  the  eigenvalues  of  H,  1 1  and  X2'  read: 


*1,2  "  Jzz«  ?  lJzz' 


(3.16) 


They  are  always  real.  Numerical  computation  is  needed  to  evaluate  them  (see  Hill. 4  and 
III.S.4).  At  t  «  0,  VP  f  they  are  both  strictly  positive,  hence  T  »  0  is  a  locus 
of  relative  minima  for  J.. 


a 


III. 3.  Tha  minimisation  algorithm 

The  algorithm  introduced  in  tlll.1  is  Implemented  on  the  discretized  counterpart  of 
the  original  minimization  problem.  The  rule  which  yields  zlc+1  e  Za<J  at  step  n.4 


SWiWi 


- ) r_4Q, 


vacuum  permittivity  e0  and  magnetic  permeability  uQ  take  on  unit  values, 

*  values  are  normalized  to  the  free  space  Impedance  zQ  :■  (Uq/Cq)1^2. 

Data  are  supplied  by  the  user  in  MKSA  units,  then  automatically  scaled.  The  computed 
functions  of  the  fields  are  printed  out  in  MKSA  units,  z  values  are  not  scaled  back. 

The  reference  frame  can  be  either  centered  at  the  source  (S)  or  at  the  mirror  image  (0) 
w.r.  to  the  impedance  plane.  Coordinates  can  be  supplied  either  in  cartesian  or  polar 
form. 

The  problem  types  which  can  be  solved  are! 

-  "direct",  i.e.  computation  of  the  discretized  counterpart  of  Jm  in  a  given 
(PIT)  domain  contained  in  Zadi 

-  "stealth",  i.e.  computation  of  z  defined  by  (3.7) i 

-  some  utility  functions:  a)  computation  of  Hessian  eigenvalues  according  to  (3.16) 
for  given  observation  point  and  z  values:  b)  computation  of  the  reciprocal 
condition  number  p  :»  Ij/Xj. 

In  connection  with  point  observation,  either  J1  (see  3.3)  or  J3  of  (3.5)  can  be 
chosen.  If  J2  (see  3.4)  is  to  be  computed,  the  user  must  specify  the  Integration 
volume.  Up  to  now  the  choice  is  restricted  to  volumes  enclosed  between  two  spherical  caps 
centered  at  0,  which  are  completely  defined  by  6  parameters,  i.e.,  see  Figure  4: 

V  {(r,9,»)|r,  «  r  *  r2i  9,  <  0  <  92i  <  *  <  e2}  .  (3.21) 

A  similar  constraint  holds  for  the  Integration  surface  E,  on  which  J4  of  (3.6)  is 
defined. 

As  soon  as  coordinates  are  read  in,  it  is  checked  that  d,  the  distance  from  the 
source  3  to  the  impedance  plane,  is  larger  than  1001  and  that  every  observation  point 
lies  in  D. 

If  a  "direct"  problem  is  set,  the  corresponding  value  of  Jm  is  computed  for  as  many 
times  as  required  by  the  user,  who  enters  the  upper  and  lower  values  of  p  and  t  and 
the  number  of  parts  into  which  both  intervals  must  be  divided. 


-19' 


If  a  "stealth"  problem  la  specified,  the  admissible  set,  (3.1)  or  (3.2)  must  be 
chosen;  then  z°  and  some  more  parameters  must  be  given.  For  simplicity,  mil. ‘miration 
in  7.^  and  in  7,^  will  be  referred  to  hereinafter  by  “unconstrained"  and  "constrained" 
respectively. 

Numerical  integration  implied  by  (3.4)  and  (3.6)  is  performed  by  the  KOROBO  routine, 
based  on  optimal  partitioning  of  the  Integration  domain,  due  to  Zakrzewska  et  al.  (14). 

Minimization  is  carried  out  as  described  in  {ill. 3.  Actually  the  rule  based  on 

k-tl 

( 3. IB)  is  chosen  in  a  fast  search  phase.  If  the  updated  value  z  of  (3.17)  is  not 
in  Za<J,  or  if  J(z  )  >  J  ,  then  a  slow  search  is  performed,  by  which  J  is  locally 
approximated  by  a  parabola,  as  described  by  Lemarechal  (15).  Function  gradients  needed  at 
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,1 


3 


every  minimization  step  are  evaluated  numerically  from  the  field  component  derivatives 
M.t.  to  z,  which  are  obtained  from  (2.12)  to  (2.19)  and  lead  to  the  procedure  shown  in 
Table  I.  Second  order  z  derivatives  of  which  we  do  not  show  here  for  simplicity 

are  also  needed  to  evaluate  ^1,2  (3*16). 

Table  1.  The  computation  of  a  field  component  z-derivatlve  as  implemented  in  the  code 


Define: 


»-  et’'r  i  ik 


Oz  \  /  *  It 


oz  n  iz 


>?0z/  Viz 


Consider  the  1st  component, 


E(jz  "  v0zCO80  '  wozsin08ln* 

«1z  *  <v1z  -  v0z  *  -  <w1z  -  w0z  *  - 

-  (vOj,8sin0  ♦  w0zecos9sin»)  •  r"1  -  ("Jocose)  *  r"1 


where  (see  (2.13)) 


UR  -  TPXjj)  •  r  -  vQ)l  •  V"  » 


wQl,  -  PYE  •  r”1  •  sinfl  1 


,  .  K-20V-1  +  Dz  -  C)  •  V"1  +  Cz)  •  (1/r2V)  . 


Explicitly: 


c  -  -O90R  -  2ReU0  ♦  2UR  -  P(T09X  +  2T0Xg  +  XggT)  +  2X(T*P  -  S"P?  )  ♦ 

♦  (R  -  P(TJX  +  XgT’  )]COS0  I 

D  -  209IUgR  +  URg  P(XTg  +  XgT)l  -  (UR  +  PXT )  •  (1  +  2UgV  1  '  Uqq  )  , 


from  which  C  ,D  are  obtained, 
z  r 


Note:  The  subscript  denotes  the  partial  derivative  w.r.t.  the  corresponding  variable: 
the  prime  (•)  division  by  sin0,  the  double  prime  (")  division  by  sin^0 . 

No  quantity  is  singular  at  9-0. 


,v'-\ 


Computational  results  are  written  out,  when  applicable,  as  printer  graphs  or  as  data 


files,  which  can  be  further  processed  by  a  graphics  package,  as  we  shall  show  in  the  next 
Section. 

III. 5.  Some  computational  results 

Remark :  All  numerical  values  are  given  in  MXSA  units  and  all  computations  have  been 
carried  out  in  single  precision  arithmetics. 

111. 5.1.  Symmetry  properties  of  J1 

The  following  symmetry  properties  of  J,  are  easily  verified  and  numerically 
checked : 

J,(9  -  0;  p  -  (1,0,0))  -  J,(9  -  0>  p  -  (0,1,0)),  VX,d,r,*,  VZ  c  Zad  ,  (3.22) 

.1,(9  -  1>  <p  -  K,  p  -  (0,1,0))  -  J,  (9  -  1;  9  -  */2>  p  -  (1,0,0)),  lA,d,r,  Vz  £  Za<J  .  (3.23) 
Polarization  effects  are  easily  evaluated  numerically.  As  expected,  if  9  yt  0,  the 
direction  of  p  affects  J,.  E.g.  for  X  «  1,  d  «  100,  r  -  200,  9  »  1  rad,  p  •  0, 
p  «  x  -  300  ohm,  p  -  (1,0,0)  yields  J,  «  1.1  IE  2  J/m3,  whereas  p  -  (0,1,0)  yields - 
J,  -  1.35E  v  2. 

111. 5. 2.  Direct  problem  solutions 

Figure  5  shows  that  the  minimum  of  J,  at  T  «  0  shifts  to  lower  values  of  p  as 
9  increases  from  0  to  1  rad. 

Table  II  gives  the  values  of  z  (ohm)  which  yield  the  minimum  and  maximum  of  J,  in 
the  specified  domain.  The  uncertainty  of  the  values  depends  on  z  domain 
discretization.  More  accurate  figures  are  obtained  by  minimizing  J,,  as  discussed  in 
till. S. 3. 


The  other  objective  functions  are  studied  in  a  similar  way.  Results  are  shown  by 
Figures  6,  7,  8. 


Table  II.  Some  results  yielded  by  the  direct  problem  solver  applied  to  W 


X  =  lm;  d  =  100mi  r  =  200m)  9=0  rad. 


10"J  4  p  <  100;  0  <  T  4  300  ohm 


M1min  at 
(J/m3)  (ohn 


Imax  at  <P 
(J/m3)  (ohm) 


0 

< 

5E  -  7 

376 

+ 

9 

0.0 

|  . 5 IE  -  3 

100. 

1/4 

< 

5E  -  7 

366 

+ 

9 

0.0 

. 47E  -  3 

100. 

1/2 

< 

5E  -  7 

328 

± 

9 

0.0 

. 38E  -  3 

100. 

3/4 

< 

5E  -  7 

270 

± 

9 

0.0 

•27E  -  3 

100. 

1 

< 

5E  -  7 

204 

± 

9 

0.0 

. 1 6E  -  3 

100. 

The  typical  CPU  time  needed  to  compute  484  values  of  or  J^  is  0.45  s.  It 

becomes  about  30  times  longer  when  J2  or  are  processed,  because  the  above  mentioned 

integration  routine  must  be  called. 

III. 5. 3.  Minimization  of  the  objective  functions 

Tables  III  and  IV  refer  to  the  minimization  of  J1.  The  terms  "unconstrained”  and 
"constrained"  have  been  defined  in  §111.4.  We  notice  that  the  machine  precision  limit  is 
achieved  in  the  "unconstrained"  case,  typically  after  30  iterations.  Given  z°  =  100  ohm, 
the  convergence  of  {zk}  towards  *exact  is  s.t.  | z*  -  *exact I /( zexact>  <  10_4  for 
k  >  16. 

In  order  to  determine  z  in  the  "constrained"  case  (p  =  t),  the  example  of  §1.4  is 

again  of  some  help.  If  we  set  a  =  n/2  in  (1.18),  (1.19)  and  consider  the  minimization 

* 

Of  •  Br,  to  which  of  (3.3)  is  proportional,  we  get: 

P  =  Zg//o"  .  (3.24) 


In  passing  we  note  that  6g/X  =  1  /ir  / 2  in  this  case.  The  z  value  just  given  must  be 
compared  to  the  computed  ones  in  line  1,  columns  3  and  5  of  Table  III  and  with  the  whole 
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The  typical  CTO  time  required  by  a  30  iteration  run  ia  0.25  s.  The  minimization  of 


J3  is  carried  out  in  a  similar  way  and  requires  comparable  CTO  times.  Some  results  are 
shown  by  Table  V.  Given  e.q.  z°  ■  1000  ohm,  the  convergence  speed  of  { Z*)  is  slightly 
higher  in  the  "unconstrained"  case  and  lower  in  the  "constrained"  one.  The  values  of  z 
computed  by  "constrained"  minimization  of  J1  and  J3  agree  to  within  3  decimal  places 
(compare  line  1,  column  8  of  Tables  IV  and  V),  which  ia  an  accuracy  test  for  the  computer 
code  and  its  implementation . 

J 2  and  J4  are  dealt  with  by  the  same  procedure,  although  for  sake  of  brevity  the 
results  are  not  presented. 


Table  V.  Minimization  of  J3  at  source  for  different  values  of  d 
X-1»jr-2dj@«»-0)P«  (1,0,0)  Vs 


computed  optimal  values  plt,Tlt  are  given  in  units  of  / u Q/e Q 


unconstrained 
z°  -  (2.652*  0.0) 


constrained 
Z°  -  (.265*  -  2*  0.0) 


d 

P* 

T* 

last 

q 

8*  -  T* 

last 

(m) 

(M/m2) 

(n.u.) 

(n.u. ) 

(W/m2) 

it.(k) 

(M/m2) 

(n.u. ) 

(M/m2) 

it.(k) 

100 

.49*  + 

11 

1.000 

.19*  -  7 

.29*  -  1 

20 

.23*  + 

12 

. 70746 

.40*  + 

11 

50 

200 

.12*  + 

11 

1.000 

.16*  -  7 

.72*  -  2 

20 

.59*  + 

11 

.71075 

.10*  + 

11 

30 

400 

.31*  + 

10 

1.000 

.158-7 

.18*  -  2 

20 

.15*  ♦ 

11 

.71082 

.25*  + 

09 

30 

800 

.77*  + 

09 

1.000 

.14*  -  7 

.45*  -  3 

20 

.36*  ♦ 

10 

.71085 

.63*  + 

09 

30 

III. 5. 4.  The  attraction  domain  of  the  minimization  algorithm 

Although  the  condition  number  criterion  (3.20)  applies  wherever  >  0,  the 
algorithm  is  shown  to  yield  a  convergent  sequence  even  if  the  starting  value  is  well 
outside  the  eubdomain  G  of  (3.19).  Minimisation  of  J.  at  source  is  successful  in  the 
dome in « 


•  1  U*  *  -  *  -  .*  .  rn* 
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